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Elasticity  of  Leucite  through  High-Temperature  Phase 

Transitions 

Donald  Isaak^,  Andy  Shen^,  Orson  Anderson',  John  Carnes' 
'University  of  California  at  Los  Angeles 
^University  of  Cambridge,  UK 
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•  Elasticity  of  Leucite 

-  Phase  transition 

.  Rutile 

-  Info  on  K '  from  high-T  RUS 

•  RUS  and  pressure  -  SiOz 
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Leucite 


Lucite 


^  White,  gray  mineral 
KAISisOe 
Important 

constituent  in  alkali 
rocks  (esp.  lava) 


TR-2 


Trademark  name  for 
a  transparent 
synthetic  resin: 
polymerized  methyl 
methacrylate 
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Phase  Transitions  in  Leucite 
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-  Leucite  - 

Variation  of  Frequencies  with  j 
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Tetragonal  Cubic 


C44—  Css—  Css—  C44 


C44°  is  the  normal  C44  without  influence  from  transition 


(See  Carpenter  et  al.,  Eur.  J.  Mineral,  10,  621-691, 1998) 


Our  Observations  with  Leucite: 

•  Wide  range  of  T  over  which  Tc  influences  C44  in  leucite 

(Compare  to  Migliori  et  al.  1993) 

•  C44  does  not  equal  C44'’  unless  d  C^VdT  >  0  (which  seems  implausible) 
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Rutile  Resonance  Spectra  (5  Modes) 
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Frequency  (Hz) 


Rutile  Resonance  Spectra  (5  Modes) 
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Frequency  (Hz) 


Temperature  Variation  of  Modes 
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Temperature  (K) 


Systematics  in  Rutile  Structure 
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Ks(M/p) 


Small  Sphere 
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Dependence  of  (dG  IdP  )j  on  Molecular  Mass  of  Gas 
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Molecular  Mass  (g) 


Basic  Equations 


B  Frequency,  ft  depends  on  G and  r 


B  Include  pressure,  Pf  and  do  algebra 
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Resonant  Ultrasound  Spectroscopy  of 
Metal-Hydrogen  Systems 
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A.  Background 

1.  Metal-hydrogen  systems 

2.  Hydrogen  motion  and  ultrasonic  attenuation 

B.  Measurements  on  a  Laves  phase  material,  TaV2,  containing 
hydrogen(deuterium) 

C.  Measurements  on  a  Ti-Zr-Ni  quasicrystal  and  a  Ti-Zr-Ni 
crystalline  approximant  containing  hydrogen 


D.  Summary 


METAL-HYDROGEN  SYSTEMS 

Many  metals  absorb  large  amounts  of  hydrogen. 

e.2.  transition  metals,  rare  earths,  actinides 


H  typically  sits  on  interstitial  sites  of  host  metal  lattice. 


Octahedral  and  tetrahedral  sites  in  the  fee  lattice 

Site  occupancy  may  be  random  (solid  solution),  or  ordered 
(hydrides). 

H  is  often  highly  mobile  and  diffuses  by  hopping  between 
interstitial  sites. 

Many  physical  effects  and  practical  applications  are 
associated  with  hydrogen  in  metals,  e.g. 

•  H>  drogen  storage 

•  Metal  h>  dride  batteries 

•  Sv»  itchable  optical  and  magnetic  properties 

•  Hx  drosen  embrittlement 


HYDROGEN  MOTION 


•  H  motion  is  studied  by  a  variety  of  techniques. 


1  -spicaJ  range  of  motion  sensed 
by  different  techniques 

NMR 

- - -  QENS 

Nfcchamcal  spectroscopy  _ 


10^  10^  108  IQio 

Characteristic  Frequency  (Hz) 


•  Ultrasonic  attenuation  arises  due  to  a  stress-induced 
rearrangement  of  H  on  the  interstitial  sites  (relaxation). 


•  Defect  symmetry  must  be  lower  than  crystal  symmetry 
for  relaxation  attenuation  to  occur  (Anelastic  relaxation 
in  crystalline  solids.  Nowick  and  Berry). 
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fflustration  of  relaxation  effects  for  sites 

on  the  face  of  a  simple  cubic  lattice 

1.  With  no  stress,  all  interstitial  sites  energetically  equal 


^  Nfetal  Atom 

O  Vacant  Interstitial  Site 
^  Occipied  Interstitial  Site 


n.  With  stress,  sites  no  longer  energetically  equal, 
rearrangement  of  interstitial  at(»ns  occurs. 
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The  loss  is  given  by 


1  AC  coTj^ 

Q  C  CO  T ^ 


.  Usually 


AC  1 

- oc  — 

C  T 


Hydrogen  can  “hop”  between  sites  by  several  different 
processes. 


Classical  Arrhenius  behavior  may  be  observed  at  high 
temperatures,  process  1. 

^R  =  '^Rc  exp(^a  /  kT) 

A  distribution  of  is  frequently  required. 

Deviation  from  classical  behavior  is  common,  especially  at 
lower  temperatures. 


I.  Sample-Transducer 
Arrangement 


11.  Spectrometer 


Transmitter 


1 

Receiver 

III.  Typical  Modes 
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Many  C-15  materials, 
including  TaV2,  absorb 
large  amounts  of  hydrogen. 

The  hydrogen  is  highly 
mobile. 

Large  number  of  tetrahedral 
sites. 


Previous  neutron  scattering  and  NMR  work  suggests  two  types 
of  motion  in  these  materials:  I.  “Slow”,  Long  range;  IL  “Fast”, 
Local 


2-D  representation  of  H(D)  in  TaV2 

Objective  of  present  study  -  Investigate  this  motion  in  a 
different  frequency  range  using  RLFS 
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Temperature  (K) 


■  1 


Can  absorb  H  up  to  H/M  =1.6 

Little  previous  study  of  H  motion  in  such  materials. 

Offers  opportunity  to  study  diffusion  in  a  non-periodic  lattice. 


•  Materials  studied 


i-phase  quasiciysial  1  i4i.5Zr4,  .Ni.t 
Loaded  with  hydrogen  to  H/M  =  1)  79 

W-phase  1/1  approximaiit  Ti44Zr4„T‘Nitf, 
Loaded  with  hA  drogen  to  H  7v1  =  0  2t  j 


H  enters  in  solid  solution. 
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Ultrasonic  loss  data  versus  temperature  for 
Ti-Zr-Ni  phase  quasicrystal. 


Temperature  (K) 


versus 


Temperature  (K) 


SUMMARY 


Evidence  for  two  types  of  motion  in  TaV2Dx  with  different 
characteristic  frequencies. 

First  ultrasonic  measurements  of  H  motion  in  quasiciystals. 
Evidence  for  a  distribution  of  activation  energies. 

Resonant  ultrasound  spectroscopy  is  a  sensitive  probe  of 
internal  motion  in  solids  in  a  frequency  range  which  is 
complementary  to  NMR  and  neutron  scattering. 
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TaV2H(D)x  Fitting  Parameters 


•  TaV.H, 

^ R  ~  ^ \  ^2 

z-j  =  r,.„  exp(£,.  /  kj) 


Eai  (eV) 

T„,(10-'W) 

(eV) 

To,(10-*sec) 

TaV2Ho34 

0.25 

3.1 

0.082 

3.6 

TaV^Hoss 

0.23 

4.8 

0.091 

1.9 

•  TaV2DQj7 

t:r  =  ^/(„exp(£,  /  kj) 

X Ro  =  3.4x10“'^  sec  =  0.25eV 
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Industrial  Applications 
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Existing  Methods 
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2  Good.  1  cracked  cu 
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More  test  types 
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Successful  Applications 
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Could  detect  chips  down  to  ~1  mm 
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Successful  Applications 
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Good  parts  often  contain  defects 


Next  Generation  Systeni(s) 
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The  Future 
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Thin  Film  Characterization 
Using  Resonant  Ultrasound  Spectroscopy 

J.  D.  Maynard,  Jason  White,  and  Jin  Hyun  So 
The  Pennsylvania  State  University 
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Micro-  (nano-)  electronics 


substrate 
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friction,  and  damping 


Use  Resonant  Ultrasound  Spectroscopy  to  measure  elastic  constants  and 
damping  for  the  film,  assuming  known  properties  of  the  substrate. 
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Measurements  are  made  for  the  substrate  with  and  without  the  film. 
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For  anisotropic  films,  the  lateral  dimensions  should  be  small  as  well 
—  Necessary  to  excite  thickness  modes  of  film 
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Electrical  leads 
to  active  area 


Beiyllium-copper 
spring  mount 
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Reproducibility  in  Mounting  and  Remounting  Sample 
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Modifications  for  Substrate  with  Film 


The  matrix  is  bloc 
four  blocks. 


PROGRAM  XYZHRQ  SUBROUTINE  FREQS  (JCVAL.  A,  YHOD,  DYDR,MA|  SUBROUTINE  XYZBLK  ( A,  HA.  WSAV) 
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GAMMA ( JG , IG ) =GAMMA { IG , JG ) 

IFd.EQ.  J)  E(IG,  JG)=F(LS,MS,NS)+DRHO*G(LS,MS,NS) 
E(JG, IG)=E(IG, JG) 


IMPLICIT  INTEGER  (I-N) 
IMPLICIT  REAL* 8  (A-H,0-Z) 
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G= ( 4 . DO /DFLOAT ( IPl * IQl * IRl ) ) * ( 1 . DO- ( 1 . DO-ALPHA) *  *  ( IRl )  ) 
RETURN 


Application:  Mats  of  Carbon  Nanotubes 
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CARBON  NANOTUBES 
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Elasticity  of  steel  and  silica  glass  spheres  under  gas  pressure 
by  RUS 

1.  Ohno\  M.  Abe\  M.  Kimura^,  Y.  Hanayama^,  H.  Oda"^  and  I.  Suzuki^ 


1  Department  of  Earth  Sciences,  Ehime  University,  Matsuyama 
790-8577,  Japan 

2  Department  of  Material  Science  and  Engineering,  Ehime  University, 
Matsuyama  790-8577,  Japan 

3  Department  of  Mechanical  Engineering,  Ehime  University,  Matsuyama 
790-8577,  Japan 

4  Department  of  Earth  Sciences,  Okayama  University,  Okayama 
700-0082,  Japan 


ABSTRACT 

Elasticity  of  minerals  under  high  pressure  are  of  great 
importance  from  geophysical  view  point.  In  order  to  develop  a 
new  method  to  measure  pressure  dependence  of  elastic  constants, 
resonant  ultrasound  spectroscopy  (RUS)  were  applied  to  the 
condition  of  gas  pressure.  We  used  a  three-layered  spherical 
shell  assembly  of  spherical  sample-thin  gas  layer-spherical 
cavity  container,  which  works  to  give  a  well-defined  boundary 
condition  in  the  analysis.  The  samples  measured  were  spheres 
of  steel  and  silica  glass  with  diameter  of  4-5  mm.  Several 
modes,  not  only  torsional  but  also  spheroidal  modes,  were 
observed  at  least  up  to  200MPa  (2kbar)  under  helium  gas 
pressure,  and  pressure  shifts  of  frequencies  were  obtained 
definitely  in  both  samples.  The  data  of  the  slope  of  pressure 
shift  of  frequency  yield  pressure  derivatives  of  bulk  modulus 
and  shear  modulus,  dK/dP  and  dG/dP,  which  are  in  good 
agreement  with  previous  data. 
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1 .  Method 


Resonant  Ultrasound  Spectroscopy  (RUS) 
with  spherical  sample 

Two  kinds  of  data  acquisition 


1)  Continuous  Wave  (CW)  method 


2)  Fourier  Transform  (FT)  method 
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Comparison  between  spectra  acquired  by  CW  method  (top)  and  FT  method  (bottom) 
sample:  silica  glass  sphere. 
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3.  Measurement  under  Gas  Pressure 


TRANSDUCERS 


X-  CAVITY 
CONTAINER 


HELIUM  GAS 


19  mm 


PRESSURE 
VESSEL 


Sample  assembly  for  measurement  under  pressure. 

~  cavity  container:  5mm  in  inner,  and  1 7mm  in  outer  diameter 
transducer:  PZT,  shear  type,  2.5MHz 
pressure  range:  1  atm  -  ZOOMPa  (2kbar) 

Eigenfrequencies  of  the  sample-cavity(gas  layer)- 
cavity  container  system  were  measured  under  gas 
pressure.  The  cavity  container  was  used  to  make  a 
well-defined  bountary  condition.  The  pressure  shift 
of  frequency  of  torsional  mode  gives  Gq',  and  those  of 

spheroidal  mode  gives  Kq'. 
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STEEL  SPHERE 


55  MPa 
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STEEL  SPHERE 
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FREQUENCY,  MHz 


STEEL  SPHERE 


Theoretical  eigenfrequencies  of  the  three-layered  spherical  structure. 
Steel  sphere  is  4.245  mm  and  the  cavity  is  1 0  mm  in  diameter. 
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Measured 


Measured 


(9  f/  (9  P  of  torsional  mode  and  G 
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df/dPof  spheroidal  mode  and  K 
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FREQUENCY  (kHz) 


SILICA  GLASS  SPHERE 
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FREQUENCY  (MHz) 


SILICA  GLASS  SPHERE 


PRESSURE  (MPa) 


Theoretical  eigenfrequencies  of  the  three-layered  spherical  structure. 
Silica  glass  sphere  is  4.764  mm  and  the  cavity  is  5  mm  in  diameter. 
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df/dP  (kHz/MPa) 


SILICA  GLASS  SPHERE 
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df/ dP  of  torsional  mode  and  G 
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(9  f/  5  P  of  spheroidal  mode  and  K 
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Silica  Glass 
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RUS  and  Rocks 


K  R  McCall 
TJ  Ulrich 

University  of  Nevada,  Reno 
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Outline 


I.  History 

II.  Experimental  Details 

III.  Numerical  Modeling 

IV.  Experimental  Results 

V.  Future  Work 
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I.  History 


A.  Empirical  Evidence:  Macroscopically 
inhomogeneous  materials  are  elastically 
nonlinear  (R  A  Guyer,  previous  talk) 


2^.0  2970.0  298O.O  2®0.0  3000.0  3010.0  3020.0 
Frequency  (Hz) 


B.  Theoretical  Models:  The  elastic 
modulus  is  nonlinear  and  hysteretic. 


o  =  Ke 

K  =  Pe  Se^  +Of(|Ac|i£)...] 

TR-3 
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C.  Hypothesis:  Reliable  measurements 
of  linear  moduli  are  the  basis  for 
quantitative  study  of  nonlinear  moduli. 


1.  Rnck  Mechanics 

<=> 


e 

destructive  &  high  amplitude 


2.  Time  of  flight 


W 


K^c^p 


path  dependent  (vs.  average  K) 


3.  Resonant  Ultrasound  Spectroscopy  ? 
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II.  Experimental  Details 


A.  RUS  System  bv 

Dynamic  Resonance  Systems.  Inc. 


control  hardware  &  analysis  software 


B.  Analysis  Software  Assumptions 

1 .  Free  boundaries  on  sample 

2.  Homogeneous  sample 

3.  Parallelepiped  sample 

C.  Conditions  for  Optimal  Results 

4.  Distinct  resonance  peaks 

5.  C„  dependence  in  low  modes 

TR-5 
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III.  Numerical  Modeling 

(to  explore  the  implications  of  the 
experimental  conditions) 

A.  Free  boundaries  on  sample 


B.  Homogeneous  sample:  Scale  of 

inhomogeneity  )  should  be  smaller  than 
wavelength  of  resonance  (A). 


L 


L>5^ 


limiting  sample  size 
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c.  ParallsIeDiped  samples:  Rock  samples  are 
cut  on  a  diamond  saw  with  limited 
accuracy,  high  probability  of  chips, 
nonparallel  sides. 

Pfiiiiirhation  Theory.  order  evaluation  of 

induced  error  (R  A  Guyer) 

u„  =  eigenmode 
Ap  =  density  change 


1  &  2  most  likely  may  be  minor  contributor 

to  error 
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Dependence  on  C 


E.  dependence  in  low  modes:  LowQ 
limits  number  of  discernable  resonances. 
Use  aspect  ratio  to  increase 
dependence. 


A 


c 


a 


b  =  A.2a 

C 

Aspect  Ratio  =  - 

a 

(basalt  and  C44) 
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IV.  Experimental  Results 

3.5 1 
3- 
2.5- 
2- 


20 

30  40  50 

60  70  80  90 

Frequency  (kHz) 

exp 

ffit  (kHz) 

1 7.988 

17.864 

19.586 

19.649 

Basalt  Sample 

23.522 

23.493 

Aspect  Ratio  =  4 

39.786 

39.692 

40.143 

40.161 

%  error  =  0.313 

41 .877 

42.071 

46.944 

46.917 

#  of  modes  =  1 0 

63.858 

63.810 

65.925 

66.142 

TR-IO 
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Basalt  Ci-,=87GPa  C44  =  32  GPa 


c 

a 

a  b  c  ^ 

(cm) 

modes 

^11  ^44 

(GPa) 

%  err 

<Q> 

1.5 

2.13  2.44  3.04  0.61 

12 

88.1  31.3 

0.632 

276 

1.7 

2.72  2.74  4.70  0.39 

10 

88.4  31.7 

1.967 

242 

2.0 

2.78  3.04  5.43  0.11 

11 

87.3  31.6 

0.709 

256 

4.0 

1.73  1.98  6.69  0.11 

10 

87.2  31.5 

0.313 

336 

Sierrj 

5  White  Granite  C 

^.|.j  =  43  GPa  ^44 

=  19  GPa 

c 

a  b  c  ^ 

modes 

P 

o 

%  err 

<Q> 

a 

(cm) 

(GPa) 

1.5 

3.68  4.16  5.64  0.66 

11 

37.6  18.7 

1.238 

146 

1.6 

2.77  3.11  4.51  0.66 

10 

43.2  18.6 

2.023 

140 

3.3 

0.96  1.20  3.17  0.44 

11 

20.8  21.3 

17.29 

149 

4.0 

1.69  2.02  6.82  0.28 

12 

37.5  19.1 

0.630 

142 
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V.  Future  Work 


A.  Same  analysis  for  cylindrical 
samples. 

(Izabela  Santos,  Senior  Thesis) 


B.  Temperature  and  Saturation  Control. 


(TJ  Ulrich,  PhD  Thesis) 


C.  Provide  baseline  for  nonlinear 
analysis. 
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Global  Symmetry,  the  Phonon 
Density  of  States,  and  RUS 

Tim  Darling  and  Albert  Migliori 

LANL 


RUS  measures  a  tiny  piece  of  the  phonon 
DOS  and  can  determine  a  global  elastic 
symmetry. 

Many  physical  properties  depend  on  the 
entire  DOS  and  the  meso-  or  microscopic 
symmetry. 

How  good  is  the  connection? 


For  when  you  are  put  into  the  Vortex  you  are  given  just 
one  momentary  glimpse  of  the  entire  unimaginable  Infinity 
of  creation,  and  somewhere  in  it  a  tiny  little  marker,  a  micro¬ 
scopic  dot  on  a  microscopic  dot,  which  says  ‘You  are  here.* 


TR-l 


149 


Harmonic  approximation: 


Replace  3N  degrees  of  freedom  by  3N 
harmonic  oscillators  (normal  modes) 

E  =  (n+1/2)hf,  n  boson  statistics 

Born-Von  Karman  (cyclic)  boundary 
conditions  in  an  infinite  solid  —  “travelling 
wave”  solutions 

No  dissipative  processes 

No  thermal  expansion 
No  temperature  dependence  of  moduli 
Isothermal  and  adiabatic  moduli  equal 

Debve  Model 

Continuum  (isotropic) 

Frequency  cutoff  by  requiring  3N  modes 
Linear  dispersion  in  all  branches 
DOS  a  f  (and  T^) 
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PLOT  OF  t"'  -  ZENER,  1948,  LOW  TEMPERATURES 

RUS  HAS  A  TRICKY  TIME  HERE,  BECAUSE 

CO  IS  GENERALLY  HIGH  (»  kHz) 

WE  SCAN  (0  OVER  A  WIDE  RANGE 

SEPARATING  OUT  CONTRIBUTIONS  FROM  DISTINCT 
MODES  IS  TRICKY  BECAUSE  OF  THE  co  VARIATION 


BUT,  WE  ARE  IN  RANGE  OF  CERTAIN  MOBILE 
INTERSTITIALS  (H) ,  AND  DOMAIN  WALLS  AND  TWIN 
AND  PHASE  BOUNDARIES  NEAR  PHASE  TRANSITIONS 
WHERE  THE  MOBILITY  IS  HIGH 
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I  Ta  polycrystal  data 
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Density  of  States  [l/THzJ 


C{T)  =  \g{w)C{a,T)da 

0 


Frequency  [T^z] 

RUS  DOS  -  Relevant  for  T  ~  few  |aK  , 
otherwise,  every  material  we  measure, 
looks  like  a  Debye  solid, 

g(a>)  a  0)^ 

-  Don’t  see  optical  branches 

-  Statistics  too  poor  to  see  anomalies  or 
contributing  acoustic  branches 
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111  zc  110  100 


DISPERSION  CURVES  FOR  Si,  Ni 


REOOCEQ  WAVE  VECTOR  COOROINATE  (5) 
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AuZn  single  crystal  =  resonance  frequenci 
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“Magic  numbers”  for  RUS  fits 


Number  of  lines  <  np*8 

(80  or  96) 

(could  be  much  lower 
depending  on  Cy  and  dim.’s) 

Number  of  lines  >  nc*5 

(10,15,25,30,45) 

rms%  <  0.5%  AMD  all  Acy  <  2% 


rms%  changes  of  <  20%  may 
not  be  meaningful 
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Contactless  Mode-Selective  Resonance  Ultrasound 
Spectroscopy:  Electromagnetic  Acoustic  Resonance 


HmoTSUGU  Ogi‘’^  and  Hassel  Ledbetter^ 
‘Osaka  University 

^National  Institute  Of  Standards  And  Technology 
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Fig.1  Typical  EMAR-measurement  setup 
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(a)OD  group 


{b)OY  group 


Lorentz  force 


(c)OX  group 


Fig. 2  Selection  of  vibration  group  with  EMAR. 
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Amplitude  (a.ii.)  Amplitude  (a.u,; 


0.3  0.4  0.5  0.6  0.7 

Frequency  (MHz) 


I 

I 


Fig.  5  Resonance  spectra  for  the  copper  monocrystal  measured  by  RUS  and  EMAR. 
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Fig.  6  EMAR  and  RUS  spe 
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icy  (MHz) 


for  the  copper  monocrystal 


Table  I  Elastic  constants  and  internal  friction  of  a  copper  monocrystal  determined  by 
EMAR  and  RUS. 


C,(GPa) 

Qii" 

(10-^) 

EMAR 

RUS 

Ref.  1 

Ref.  2 

EMAR 

RUS 

C„ 

167 

168.7 

169.7 

168.4 

1.14 

1.74 

C|2 

120.9 

121.6 

122.6 

121.4 

0.3 

1.34 

C44 

74.64 

75.46 

74.49 

75.4 

2.59 

2.46 

B=(C„+2C,2)/3 

136.2 

137.3 

137.3 

137.1 

0.64 

1.5 

C-(C, -C,j)/2 

23.05 

23.52 

23.35 

23.5 

3.34 

2.8 

218.6 

220.6 

220.6 

220.3 

1.4 

1.87 

^1  1  *^2  2+4044)73 

235.7 

237.9 

237.6 

237.7 

1.47 

1.9 

Cni.arb“(Gn“Cl2“^^44)/3 

40.25 

40.84 

40,54 

40.8 

2.88 

2.59 

Ct 

199.6 

201.6 

- 

- 

1.34 

1.84 

G 

47.54 

48.22 

- 

- 

2.83 

2.57 

E 

127.8 

129.5 

- 

- 

2.6 

2.46 

V 

0.3437 

0.3428 

- 

- 

- 

1.  H.  Ledbetter  and  S.  Kim,  Personal  communication,  unpublished. 

2.  R.  Hearmon,  Adv.  Phys.,  5,  323  (1956). 
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Table  II  Elastic  constants  (GPa)  and  internal  friction  of  SiCf/Ti-6AI-4V  crossply 
composite  measured  by  EMAR,  RUS,  and  pulse-echo  methods. 


EMAR 

RUS 

pulse  echo 

EMAR 

EMAR* 

Q-'ydO-^) 

Cn 

225.4 

230.4 

- 

221.7 

222.1 

7.79 

C33 

191.6 

185.8 

188 

182.6 

183.3 

15.6 

C44 

54.58 

54.98 

53.7 

54.6 

54.63 

9.33 

Cgg 

51.00 

51.73 

- 

53.02 

52.08 

6.33 

C12 

80.82 

85.14 

- 

74.2 

74.36 

16.1 

^13 

79.16 

80.9 

- 

71.29 

71.65 

23.7 

En 

180.6 

182.3 

182.7 

183.1 

3.54 

E33 

149.7 

142.1 

146.8 

147.3 

9.59 

B 

124.5 

126.7 

- 

117.7 

118.1 

15.2 

*  Determined  from  OD  and  OY  vibration  groups. 
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Fig.  9  EMAR  measurement  setup  at  elevated  temperatures. 
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Relative  Amplitude 
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Fig.  10  EMAR  spectra  of  OY  group  at  various  temperatures  for 
SiC,/Ti-6AI-4V  crossply  composite. 
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33,  B  (GPa)  C44,  Q/,  (GPa)  C,i,  C33  (G 


300  400  500  600  700  800  900  1000 

Temperature  (K) 


Fig.  11  Temperature  dependence  of  elastic  constants  of  SiCf/Ti-6AI-4V 
crossply  composite. 


Temperature  (°C) 


Fig.  12  Temperature  dependence  of  internal  friction  of  a  OY  mode  of 
SiC(/Ti-6AI-4V  crossply  composite.  Arrows  indicate  internal  friction  after 
cooling  down  to  room  temperature. 
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Elastic  Coefficients  and  Internal  Friction 
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Si  elastic-stiffness  coefficients 
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Source:  S.  Kim,  H,  Ledbetter,  NIST — Boulder.  26  MAY  1999 
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From  the  Koehler-Granato-Lucke  stretched-vibrating- 
string  model,  we  have  for  low  frequencies  (in  modified 
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TR-11 


Fig.  14-1.  Internal  friction  as  a  function  of  strain  amplitude  for  an  undeformed 
copper  crystal  (curve  A)  and  after  application  of  stresses  of  60,  120,  and  150  psi,  respec¬ 
tively  (curves  B,  C,  and  D).  Measurements  in  longitudinal  vibration  at  about  30  kHz. 

(From  Read,  1941.) 
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Elastic  Constants  and  Internal  Friction  of  Quartz 

and  Langasite 
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*Now  at  Lucent,  Denver 
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Langasite  Quartz 


Phase  transform. 

573°C 

Q-i  (10-5)** 

17.3 

12.6 

^11 

168 

39 

€33 

460 

41 

en  (C/m^  ) 

-0.43 

-0.173 

ei4  (C/m^) 

0.11 

0.04 

4.2 

2.5 

*Equivalent  designations:  P321,  D3 ,  hP9. 
**First  torsional  mode  of  a  cylinder. 
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Table  1.  Langasite  and  Quartz:  basic  properties 


Langasite 

Quartz 

Formula 

La3Ga5SiOi4 

SiOj 

Mass  density  (g/cm^) 

5.744 

2.6500 

Point  (space)  group 

Trigonal  32  (150)* 

Trigonal  32  (150) 

Atoms/cell 

23 

9 

C„  (GPa) 

186.1 

87.05 

C33  (GPa) 

259.8 

106.9 

Vn  (cm/ps) 

0.5692 

0.5731 

V33  (cm/|is) 

0.6725 

0.6328 
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Quartz  elastic-stiffness  coefficients 
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Quartz  and  Langasite  elastic-stiffness  coefficients 


Quartz 

Langasite 

Cjj  (GPa) 

P 

1 

Cii  (GPa) 

p  (g/cm^) 

2.650 

5.744 

^11 

87.10 

-0.09 

186.1 

0.18 

C33 

106.1 

0.66 

259.8 

0.27 

C44 

58.26 

-0.07 

53.87 

0.07 

Cee 

40.17 

-0.17 

42.68 

0.03 

C12 

6.76 

100.7 

Ci3 

12.17 

93.50 

^14 

-18.07 

14.08 

rms  (%) 

0.19 

0.26 

Cl  (GPa) 

100.7 

0.11 

202.4 

0.22 

B  (GPa) 

37.80 

0.41 

133.0 

0.31 

G  (GPa) 

47.18 

-0.07 

52.00 

0.05 

E  (GPa) 

99.94 

0.07 

138.0 

0.08 

V 

0.0593 

0.3271 

21  June  1999 
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Langasite 
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Measured  Resonance  Frequencies  (kHz) 
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RUS  STUDIES  OF  CRYPTO-CLATHRATES: 
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NCPA,  University  of  Mississippi 


ORNL: 

B.C,  Sales 

D  Mandrus 

B.  C.  Chakoumakos 

KUL: 

C  Laermans 


LANL: 

A.  Migliori 
T.  Darlina 
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Filled  skutterudites:  RM4X12 

R=Sr,  Ba,  La,  Ce,  Pr,  Nd, 
Eu,Gd,  U,  Th 

M=Ru,  Fe,  Os 

X=P,  As,  Sb 

R-atoms  “rattle”  in  atomic  cage 
"  reduce 

excellent  example  of  “electron-crystal 
phonon  glass” 

New  class  of  thermoelectric 
materials 


Clathrates:  XqE^q 

X=  Na,  K,  Rb,  Cs,  Sr,  Ba 
E  *  Si,  Gg,  Sn,  Al,  Ga 


210 


(W/mK) 


CRYPTO-CLATHRATbS: 


RUS  on  filled  and  unfilled  skutterudite 
antimonides; 

CoSbj 

polycrystalline  — ►  and  C44 

T:  5  -  300  K 

US  absorption  (pulse-echo)  on  Ge-clathrates 

Sr8Gai2®®30 
single  crystal 


0^4  (10  N/m'^)  044(10  N/m 
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elastic  constants  for  two-level  system 

c=  F=  -NAkBTIn(1  +e‘^) 

A  =  Ag-l-Ae 


c  =  c  (background) 

+  c(TLS#1,  Aoi=50  K) 

+  c  (TLS  #2,  Ao2=  200  K) 

c  (background)  =  Varshni-fit  for 

unfilled  sample 
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0^4  (10"  N/m‘) 


0  100  200  300 

Temperature  (K) 


C44(model)  =  C44(background) 

+  C44(TLS#1,  Aoi=50  K) 

+  C44  (TLS  #2,  Ao2=  200  K) 

tOoXurc.  ^ 
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100  200 
Temperature  (K) 


(model)  =  (background) 

+  Cii(TLS  #1,  Aoi=  50  K] 


+  Cii  (TLS  #2,  Ao2= 


lie 


Sr.Ga.^Ge 


(Luo/gp)  uojidjosqv  sn 
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Temperature  (K) 


2  local  modes  in  fiiied  skutteruaites 


tunneling  states  in  G 


perfect  crystals  with  glasslik 
elastic  properties 


future  work: 

RUS  study  of  Ge-clathrate 
US  attenuation  of  skutterudites 


Comparison  of  radiation  and  scattering  mechanisms  for  objects 
having  Rayleigh  wave  velocities  greater  than  or  smaller  than  the 
speed  of  sound  in  the  surrounding  water 

Philip  L.  Marston 
Florian  J.  Blonigen 
Brian  T.  Hefner 
Karen  Gipson* 

Scot  F.  Morse** 

Dept,  of  Physics 
Washington  State  University 
Pullman,  WA  99164-2814 

Supported  by  the  Office  of  Naval  Research. 

*Physics  Department,  University  of  Puget  Sound,  1500  North 
Warner,  Tacoma,  WA  98416 

**Physical  Acoustics  Branch,  Code  7136,  Naval  Research 
Laboratory,  4555  Overlook  Ave.,  SW,  Washington,  DC  20375 

•  Fluid  loading  mechanisms  for  metallic  objects  in  water:  review 
and  recent  progress  (Morse,  Gipson,  Blonigen). 

•  Reflection  from  plastics,  fluid  loading,  and  the  subsonic 
Rayleigh  wave. 

•  Novel  backscattering  enhancements  for  plastic  (and  rubber) 
objects  in  water  (Blonigen,  Hefner). 
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cylinder 
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form 

function 


reflectance 
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Space-Time  Approach  to  Non-Relativistic 
Quantum  Mechanics 

R.  P.  Feyxman 

Cornell  University,  I  theca,  York 
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=  J  (18)  ^  ^  ^ 


7.  DISCUSSION  OF  THE  WAVE  EQUATION 
The  Classical  Limit 

This  completes  the  demonstration  of  the  equiv¬ 
alence  of  the  new  and  old  formulations.  We 
should  like  to  include  in  this  section  a  few  re¬ 
marks  about  the  important  equation  (18), 

This  equation  gives  the  development  of  the 
wave  function  during  a  small  time  interval.  It  is 
easily  interpreted  physically  as  the  expression  of 
Huygens’  principle  for  matter  waves.*  In  geo¬ 
metrical  optics  the  rays  in  an  inhomogeneous 
medium  satisfy  Fermat's  principle  of  least  time. 
We  may  state  Huygens*  principle  in  wave  optics 
in  this  way:  If  the  .amplitude  of  the  wave  is-' 
known  on  a  given  suriace,  the  amplitude  at  a 
near  by  point  can  be  considered  as  a  sum  of  con¬ 
tributions  ^rom  all  points  of  the  surface.  Each 
contribution  is  delayed  in  phase  by  an  amount 
proportional  to  the  time  it  would  take  the  light  to 
get  from  the  surface  to  the  point  along  the  ray  of 
least  time  of  geometrical  optics.  W*e  can  consider 
(22)  in  an  analogous  manner  starting  with 
Hamilton’s  first  principle  of  least  action  for 
classical  or  "geometrical”  mechanics.  If  the 
amplitude  of  the  wave  is  known  on  a  given 
"surface,”  in  particular  the  "surface”  consisting 
of  ail  X  at  time  t,  its  value  at  a  particular  nearby 
point  at  time  t~-t,  is  a  sum  of  contributions  from 
all  points  of  the  surface  at  t.  Each  contribution  is 
delayed  in  phase  by  an  amount  proportional  to 
the  action  it  would  require  to  get  from  the  surface 
to  the  point  along  the  path  of  least  action  of 
classical  mechanics.^®  TR-lO 
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^mm 


Observed  far-field  backscattered  form  function 


Karen  Gipson 
1998  Ph.D.  thesis  ^ 


^  1-5 
1 

0.5 


26  28  30 

Y(deg) 


«  ^Cu 


frequency=  1 .03  MHz 

'N  i  r  „  .  , 

5-  ^  e*^'' 


data 

theory 

shifted  theory 


theory  peak 

%  discrepancy 

3.45 

25% 

3.20 

5% 

620  kHz  2.60 


1.03  MHz  3.05 


Fluid  Loading  Properties  Relevant  to  Scattering: 

c  (m/s)  p  (g/cc) 

1.  Water  1490  1.0 

2.  Sandy  Sediment  1800  1.8 

(fluid  model) 


Elastic  Parameters  (Solid  or  Shell  "Target"): 


CR 

cs 

CL 

P 

3.  Stainless  Steel 

2904 

3141 

5675 

7 

4.  Solid  Polymers: 

A.  PMMA 

1250 

1340 

2690 

1.19 

Polystyrene 

1075 

1150 

2400 

1.05 

B.  RTV  Rubber 

- 

- 

1080 

1.02 

What  are  some  mechanisms  for  enhanced  backscattering 
from  TYPE  4(A)  or  (B)  targets  with  TYPE  1  or  2  loading? 
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PLASTIC  PROPERTIES 


GLASS  TRANSITION  -  a  second  order  phase  transition 
from  the  glassy  state  to  the  rubbery  state  having  a 
negligible  shear  modulus. 

GLASS  TRANSITION  TEMPERATURE  (Tg)  - 

the  temperature  at  which  the  glass  transition  occurs 
in  the  low  frequency  limit. 


Tg  (C) 

PMMA  105 

Polystyrene  1 00 

•  While  the  shear  wave  velocity  Cs  is  relatively  dispersive, 
for  T  ~  20  C  at  low  ultrasonic  frequencies, 

AC3  «  0.01  C3 

for  the  frequency  range  in  our  experiments. 

•  With  k3  =(co/c3)(1  +  i63) 

Rl  =  (co  /  C|_)(1  +  i5|_) 

5s  « 1  and  5|_  « 1  at  low  ultrasonic  frequencies  and 
do  not  depend  on  co.  Example:  PMMA  5s  ~  0.0053 
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Fig.  4.  Primary  and  secondary  Descartes  rays  projected  onto  the 
base  plane  for  a  circular  ice  cylinder  with  n  -  1.31  and  illuminated 
at  y  =  20°.  As  shown  for  this  case  in  Fig.  3,  the  Descartes  ray 
scattering  angle  for  p  =  2  exceeds  Gjt)  for  p  =  3.  The  projected 
specular  reflection  and  bisector  of  the  cylinder  are  shown  from  the 
perspective  of  a  distant  observer  at  the  of  the  primary  rainbow. 
The  meridional  plane  that  contains  the  cylinder’s  axis  and  the 
incident  wave  vector  is  shown.  Glare  points  associated  with  spec¬ 
ular  reflection  were  not  as  bright  as  those  for  the  Descartes  rays  for 
situations  in  which  the  glare  points  from  Descartes  rays  were 
visible. 


the  base  plane  for  the  cases  ofp  =  2  and  3,  which  are 
analogous  to  primary  and  secondary  rainbows.  The 
filled  circle  shown  is  the  location  of  the  corresponding 
ray-optics  discontinuity  in  Fig.  5  of  Ref.  8  for  the 
specific  case  of  y  =  45*^.  This  comparison,  as  well  as 
laboratory  experiments,^®  gives  support  to  this  anal¬ 
ysis.  As  the  tilt  angle  increases,  the  p  =  2  caustic 


evolves  in  such  a  way  that  6p  reaches  180°  at  a  crit¬ 
ical  tilt  angle  given  by  taking  n'  =  2  in  Eq.  (1). 
This  corresponds  to  a  merging  in  the  meridional  scat¬ 
tering  plane  of  the  Airy  caustics  from  both  sides  of  the 
cylinder  (see  Appendix  A). 

For  small  values  of  tilt,  the  evolution  of  the  values 
of  0x)  and  the  locations  of  the  associated  rays  are  more 
closely  related  to  the  observations  of  sunlit  icicles. 
Figure  4  shows  the  locations  of  the  Descartes  rays 
projected  onto  the  base  plane  for  p  =  2  and  3  as  well 
as  the  specular  ray  (which  may  be  assigned  the  value 
p  =  0),  which  has  the  same  scattering  angle  as  0£)  for 
the  primary  rainbow  ray.  For  the  tilt  angle  y  of  20° 
used  in  this  figure,  an  observer  positioned  to  view  the 
glare  of  the  primary  rainbow  ray  would  be  in  the 
shadow  zone  of  the  secondary  rainbow  (see  Fig.  3)  so 
that  no  rays  having  two  internal  reflections  would  be 
visible.  Conversely,  an  observer  positioned  to  view 
the  glare  of  the  secondary  rainbow  ray  would  be  in 
the  shadow  zone  of  the  primary  rainbow. 

3.  Caustics  and  Glare  Points  Created  by  Sunlit  Icicles 

Sunlit  icicles  were  observed  visually  and  photograph¬ 
ically  in  Pullman,  Wash.,  during  February  1996. 
Most  of  the  observations  and  supporting  photographs 
were  taken  in  midmoming  with  the  Sun  position  cal¬ 
culated  to  be  20°  to  25°  above  the  horizon.  With  the 
axis  of  the  icicle  approximated  as  vertical,  the  corre¬ 
sponding  values  for  y  in  the  cylinder  model  range 
from  20°  to  25°.  Figures  5(a)  and  5(b)  show  color 
photographs  with  the  camera  focused  on  the  icicle 
such  that  the  camera  is  on  opposite  sides  of  the  me- 


(a)  (b)  (c) 


Fig.  5.  (a)  Descartes  glare  point  for  the  primary  Airy  caustic  of  an  icicle  is  visible  as  the  bright  pa^  on  the  right-hand  side  of  the  icicle 
near  its  center.  The  weak  glare  point  visible  on  the  left-hand  side  is  for  a  twice-reflected  and  twice-refracted  ray.  (b)  The  same  icicle 
viewed  from  the  other  side  of  the  meridional  plane  shows  a  line  of  primary  Descartes  glare  points  on  its  left-hand  side  with  a  hne  of  specular 
glints  to  the  right  of  an  imaginary  center  line  (shown  as  the  cylinder’s  bisector  in  Fig.  4).  The  lower  of  the  two  very  bright  glare  pomts 
has  a  green  hue.  (c)  Perspective  as  in  (b)  but  with  the  camera’s  focal  plane  conjugate  to  the  far-field  scattering.  ’The  glare  pomts  are  now 
out  of  focus,  and  some  of  the  Descartes  glare  points  appear  to  be  associated  with  colored  regions. 

TR-15 
236 


20  March  1998  /  Vol.  37.  No.  9  /  APPLIED  OPTICS  1553 


ridian  formed  by  kj-  and  the  cylinder’s  axis.  Glare 
points  evident  in  these  photographs  were  also  easily 
observed  by  eye.  The  causes  of  the  glare  points  were 
identified  when  the  head  of  the  observer  moved  and 
by  comparisons  with  ray  diagrams  as  shown  in  Fig.  4. 
Note  that  Fig.  4  is  drawn  for  the  perspective  shown  in 
Fig.  5(b)  such  that  the  Descartes  ray  appears  on  the 
left  whereas  a  glint  from  external  specular  reflection 
appears  on  the  right-hand  side  of  the  bisector.  Cor¬ 
responding  features  can  be  identified  in  Fig.  5(b). 
The  lowest  bright  glare  point  has  a  green  hue. 

The  glare  points  (or  point)  from  the  two-chord  ray 
in  Fig.  5(b)  disappeared  when  6  was  reduced  below  a 
value  that  was  interpreted  to  be  6^  because  the  tran¬ 
sition  appeared  to  be  associated  with  the  Descartes 
ray.  In  support  of  this  interpretation,  the  glare  was 
observed  to  be  especiaUy  bright  at  this  Airy  transition 
when  compared  with  the  glare  of  the  specular  reflec¬ 
tion.  Furthermore,  observations  with  a  camera  defo- 
cused,  so  the  focal  plane  was  conjugate  to  the  far-field 
scattering,  manifested  the  color  sequence  of  the  rain¬ 
bow  in  the  defocused  image  of  some  two-chord  glare 
points  but  not  for  the  specularly  reflected  glint.  (This 
viewing  method  displays  the  angular  dependence  of 
the  scattered  light.  Other  examples  of  out-of-focus 
viewing  are  discussed  on  pp.  231-234  of  Ref.  7.)  Fig¬ 
ure  5(c)  is  a  far-field  image  for  the  camera  perspective 
as  in  Rg.  5(b)  and  Fig.  4.  Vertical  colored  bands  ap¬ 
parently  associated  with  the  Airy  caustic  of  the  pri¬ 
mary  rainbow  rays  are  visible  slightly  below  the  center 
of  the  photograph,  with  the  red  side  (and  the  corre¬ 
sponding  shadow  region  of  the  caustic)  appearing  on 
the  right-hand  side  of  the  bands.  It  was  confirmed 
with  a  diagram  based  on  elementary  lens  theory  that 
the  shadow  side  in  the  defocused  photograph  [Fig.  5(c)] 
agreed  with  predictions  in  which  the  illuminating  rays 
correspond  to  the  primary  rainbow  rays  of  the  icicle. 
The  limited  extent  of  the  colored  bands  in  Fig.  5(c) 
appears  to  be  the  result  of  vignetting  due  to  the  lens 
aperture.  With  the  observer  in  the  two-ray  region  of 
the  primary  rainbow,  it  was  difficult  to  resolve  the 
distinct  glints  of  the  associated  rays  (labeled  1  and  2  in 
Fig.  2);  however,  that  may  be  more  a  deficiency  of  the 
observer’s  vision  (the  author’s)  than  of  the  actual  op¬ 
tical  phenomenon. 

Although  it  was  not  feasible  to  measure  for  nat¬ 
ural  icicles,  some  information  about  relative  positions 
of  Descartes  angles  is  avsdlable  because  glare  points 
that  have  the  expected  properties  of  a  secondary  or 
p  =  3  Descartes  ray  were  occasionally  visible  by  care¬ 
ful  positioning  of  the  observer.  Such  a  glare  point  is 
visible  on  the  left-hand  side  of  the  icicle  in  Fig.  5(a). 
This  glare  point  is  slightly  below  a  primary  Descartes 
ray  glare  point  visible  on  the  right-hand  side.  Al¬ 
though  Fig.  4  shows  that  the  primar\^  and  the  sec¬ 
ondary  glare  points,  if  visible,  should  be  on  opposite 
sides  of  a  cylinder,  the  recording  of  both  glare  points 
with  a  single  camera  position  merits  some  discussion. 
This  is  because  inspection  of  Figs.  3  and  4  shows  that 
the  secondary  rainbow  lies  in  the  shadow  region  of 
the  primaiy  rainbow  of  a  cylinder  except  for  small 
values  of  7,  namely,  7  <  4.9®  for  n  =  1.31.  The 
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simultaneous  observation  of  the  primary  and  the  sec¬ 
ondary  glare  points  as  in  Fig.  5(a)  may  be  indicative 
of  deviations  from  a  circular  cross  section  or  a  cylin¬ 
drical  shape,  both  of  which  are  visible  in  Fig.  5  and 
from  direct  inspection  of  the  icicle.  It  is  noteworthy 
here  that  the  vertical  spread  of  the  glints  from  the 
p  =  2  ray  also  indicates  de\dations  from  a  cylindrical 
shape,  since  it  is  to  be  expected  that,  for  perfect  cyl¬ 
inders,  the  vertically  displaced  glare  points  would 
spread  vertically  and  merge  into  vertical  lines. 

Although  Fig.  5  and  the  observations  described  are 
based  on  an  icicle  that  has  a  diameter  of  approxi¬ 
mately  1  cm,  the  qualitative  features  of  the  primary 
Descartes  ray  were  confirmed  with  observations 
based  on  clusters  of  smaller  icicles  in  December  1996. 

4.  Discussion  and  Conclusions 

The  observations  summarized  in  Section  3  show  that 
rays  within  icicles  can  produce  bright  glare  points 
that  have  the  expected  properties  of  primaiy  and 
secondary  Descartes  (or  rainbow)  rays.  The  obser¬ 
vations  were  first  made  with  the  unaided  eye  and  are 
supplemented  by  color  photographs,  as  shown  in  Fig. 

5.  The  observations  were  made  without  the  aid  of  a 
polarizer  and  with  solar  illumination  unobscured  by 
clouds. 

A  model  was  presented  for  how  the  projected  scat¬ 
tering  angle  for  primary  and  secondary  Descartes 
rays  varies  with  the  tilt  angle  7  of  the  illumination 
relative  to  the  axis  of  the  icicle  for  the  situation  in 
which  the  icicle  has  the  shape  of  a  circular  cylinder. 
This  model  neglects  the  taper  and  the  shape  varia¬ 
tions  of  real  icicles.  For  a  sufficiently  small  value  of 
7,  corresponding  to  nearly  horizontal  illumination, 
the  model  suggests  that  a  sufficiently  distant  ob¬ 
server  could  lie  within  the  two-ray  region  of  both  the 
primaiy  and  the  secondary  caustics.  The  observa¬ 
tions  were  made  with  7  only  somewhat  larger  than 
the  crossover  value  (predicted  to  be  at  or  near  4.9®)  so 
that  the  apparent  simultaneous  observation  of  pri¬ 
mary  and  secondary  rainbow  glare  points  in  Fig.  5(a) 
may  be  due  to  deviations  from  the  simplified  shape 
assumed. 

The  above-mentioned  model  also  predicts  that  at  a 
critical  value  of  tilt  7^  (predicted  to  be  at  60.75®  for 
ice)  the  primary  rainbow  Descartes  rays  from  both 
sides  of  the  cylffider  merge  in  the  meridional  plane. 
The  expression  for  7^.,  Eq.  (A3),  may  be  derived  either 
with  the  Bravais  effective  refractive  index,  Eq.  (1),  or 
from  a  direct  calculation  of  the  appropriate  principal 
curvatiu-e  of  the  refracted  internal  wave  front  in  the 
meridional  plane,  Eq.  (A5).  Although  it  is  not  cur¬ 
rently  known  whether  the  caustic-merging  transition 
(CMT)  is  observable  for  naturally  sunlit  icicles,  if  it  is 
observable  for  icicles,  the  transition  glare  point 
should  be  very  bright.  This  consequence  of  the  en¬ 
hanced  focusing  is  discussed  in  Appendix  A.  The 
merging  and  the  general  brightness  enhancement  of 
the  two-chord  glare  point  were  both  confirmed  by  the 
viewing  of  a  long  tilted  polished  glass  rod  illuminated 
with  light  from  a  distant  slide  projector.  Unfortu¬ 
nately,  conditions  favorable  for  such  observations 
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Scattering  observations  for  tilted  transparent 
fib^:  evolution  of  Airy  caustics  with  cyRnder  tilt 
and  the  caustic  merging  transition 

Catherine  M.  Mount,  David  B.  Thiessen,  and  Philip  L.  Marston  APPLIED  OPTICS  /  Vol.  37.  No.  9  20  March  199S 
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New  Backscattering  Enhancement  Mechanisms  for  Solid  Polymer 
Targets 


1 .  For  filled  targets  (such  as  solid  cylinders):  "New"  transmitted 
rays  for  generating  flat  scattered  wavefronts. 

Examples:  Subsonic  shear  waves  (PMMA  or  Polystyrene) 
Subsonic  longitudinal  waves  (RTV) 

Geometry:  Tilted  bluntly  truncated  cylinder 

Caustic  merging  transition  at  critical  tilt: 


Yc  =  arccos 


,1/2 


.  N>1, 


Refractive  Index  N  =  c  /  Cj  or  c  /  C|_ 

Marston,  Applied  Optics  37, 1551-1556  (1998) 

Blonigen  and  Marston,  JASA  (abstract)  102,  3088  (1997) 
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w,  ^  ^  =  J  ds  exp  [+  i  (s‘‘+  w^sVw.si] 


•“OO 


l^co.6-)  =  irc-5;)e*'? 

8 e.cW scattered  Pressure  : 

P=  P.- 

cohere  'f  -  "Foi^n^  "For%c*t:ion 

hc^rcX  sphere  •  |-f|  =  1 

tLiHed  cy\5ndler  :  On-fUe  Cdas*Hc  W^=W2^“  0  ^ 

1^'=  |rC!,)lBl 
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New  Backscattering  Enhancement  Mechanisms  for  Solid  Polymer 
Targets 

2.  Tunneling  enhancements  for  “new”  subsonic  waves  .  ' 

•  Example:  Rayleigh  wave  on  PMMA  or  Polystyrene  becomes 
subsonic  with  respect  to  water  or  sediment  loading. 

•  Ray  Theory:  Same  as  Zhang,  Sun,  &  Marston  (1 992)  but  with 
Rayleigh  wave  parameters  and  attenuation. 

•  Experimental  confirmation  for  solid  spheres:  B.  T.  Hefner. 


Rayleigh  wave  velocity  for  an  acrylic  sphere 


in  water. 


l->  3.  New  coincidence  frequency  and  backwards  wave 
■  enhancements  for  sheiis. 
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Summary 


•  Quantitative  leaky  ray  theory  has  been  confirmed  for  meridional 
ray  backscattering  enhancements  for  truncated  steel  shells  and 
solid  cylinders  (Morse  and  Gipson). 

•  In  addition  to  leakage  by  supersonic  ao  waves,  reflection  from  the 
end  of  a  shell  near  a  mode  threshold  causes  radiation  due  to  the 
collective  motion  of  a  cut-off  wave  (Morse). 

•  Plastic  targets  In  water  or  sediment  require  rethinking  of 
enhancement  mechanisms. 

•  For  blunt  tilted  solid  cylinders  of  polystyrene  and  RTV  a 
backscattering  enhancement  due  to  the  merging  of  farfield 
caustics  of  transmitted  rays  has  been  observed  (Blonigen). 

•  PMMA  spheres  were  observed  to  give  a  backscattering 
enhancement  due  to  tunneling  to  a  subsonic  Rayleigh  wave 
(Hefner). 

•  We  modeled  the  resulting  Rayleigh  wave  resonance  scattering 
by  extending  previous  ray-tunneling  models  of  subsonic  ao 
waves  on  steel  shells  [Zhang,  Sun,  &  Marston  JASA  (1992)]. 
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Subsonic  Rayleigh  wave  resonances  on  solid  polymer  spheres  in 
water  and  backscattering  enhancements  associated  with 
tunneling:  experiments,  models,  and  the  relative  significance  of 

material  and  radiation  damping 


Brian  T.  Hefner 
and 

Philip  L  Marston 
Dept,  of  Physics 
Washington  State  University 
Pullman,  WA  99164-2814 


Supported  by  the  Office  of  Naval  Research. 
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Backscattering  from  solid  spheres 


Material  properties: 


Longitudinal  velocity 
(mm/|is) 

Shear  velocity 
(mm/ps) 

Rayleigh  velocity 
(mm/ps) 

Density 

(g/mm') 

Stainless  Steel 

5.675 

3.141 

2.903 

7.57 

Acrylic  (PMMA) 

2.690 

1.340 

1.250 

1.19 

Water 

1.479 

1.00 

Coupling  of  incident  field  to  subsonic  Rayleigh  wave: 


B, 


Partial  wave  series  (PWS)  solution  : 


ika 


Sommerfeld  Watson  Transform  (SWT): 


Ov,  M  =  0 

V,  =a,  +(p, 

Cl  _  X 

c  oc^  +  ^ 


Cl :  phase  velocity  for  the  I  th 
surface  wave 

Pi :  damping  coefficient 


Rayleigh  wave  contribution: 

fmR  =  ~^R  “  27tmP/j '+  i{r\ii  -mK-\-  2nmx  cj )] 

y.  -G;;exp(-7cP^+iTl;;) 

^  [l  +  exp(-27ip^ +i27D:c/c/j)] 

G^|~87ip;jC/c^  =7t/:ac/c^ 
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Incorporating  material  absorption  into  the  backscattering 

formulations. 


Assume  both  longitudinal  and  shear  wave  velocities  have  linear  absorptions 

ki  kj^(l  +  ia^ )  =  0.0034 

k,=>k,(l  +  ia,)  For  Acrylic:  „^^o  o053 

In  the  exact  Ifl- 


Ifi  for  Solid  Aaylic  Sphere 


damped  Ifl  for  Solid  Aaylic  Sphere 
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10-cycle  tone  burst  backscattered 
from  acrylic  sphere. 


Sheet  Source 


\ 

Acrylic 
Sphere 
(r  =  r) 


For  calculated  resonance 
at  1 5.9  kHz 

Measured:  Q  =  9.31 


14.5  kHz  15.0  kHz 


0  0.2  0.4  0.6  0.8  1  1.2  0  0.2  0.4  0.6  0.8  1  1.2 


16.5  kHz  17.0  kHz 
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Absorption  in  the  approximate  Rayleigh 
wave  contribution. 

G^l  ~  STCp^c/  depends  on  the  radiation  damping 
parameter,  pR 

In  the  remaining  expression  for  fmi,  the  damping  parameter 
must  include  both  radiation  damping  and  material 
absorption. 

The  equation  {x)  =  0  is  then  solved  with  the  complex 

longitudinal  and  shear  wavenumbers  to  yield  +/P'^ 

and 


oo 


m=0 


exp(-7ip'^  +/T1^) 

[1  -  exp(-27cP^  +  ilnx  cj  )] 


Possible  to  approximate  oc'^  and  P'^  by  using  results  for 
fluid-loaded  PMMA  half  space: 

+  =>  Yr  =  0.00149 


Using  this  value  and  the  values  calculated  for  the  sphere 
without  material  absorption: 


^  R{approx) 
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Rayleigh  wave  velocity  for  an  acrylic  sphere 

in  water. 


Rayleigh  wave  damping  for  an  acrylic  sphere 

in  water. 


ka 
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and  exact  form  function 
absorption 


Measurements  of  form  function 
for  acrylic  sphere. 


Transducer 


/ 

Hydrophone 


Acrylic 

sphere 


1—  68.5  cm  — 
1.87  m - 


30-cycle  tone  bursts 
from  a  3-1  composite 
transducer. 

The  response  was 
normalized  by  measuring 
the  signal  when  the 
hydrophone  was  placed 
at  the  sphere  location. 


Data  taken  for  two  spheres  ^  a  =  1 .27  cm 
in  the  frequency  range 

30  -  60  khz  ♦  a  =  2.54  cm 


Comparison  of  mesured  and  calculated  |f| 
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Calculation  of  quality  factors  for 
Rayleigh  wave  resonances  of  a 
fluid-loaded  PMMA  sphere. 

Quality  factors  of  resonances  can  be 
calculated  from  decay  rate  of  echo 
amplitudes  (L.  R.  Dragonette,  et.  al.  1996) 

mthecho  ^  A^-B  e\^{-TPRCgRlc) 

T=(ct-  R)/a 
Q  «  ka/(2pR  CgR/c) 
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Target  Strength  for  PMMA  spheres. 

Thin  fluid-filled  shells 
typically  used  as 
passive  sonar  targets. 

Glory  scattering 

Cw<Cf 

CFC's  have  typically  been  used  as  filler  fluid  however 
alternatives  are  now  required  due  ban  on  CFC's. 
(G.  Kaduchak  and  C.  M.  Loeffler,  1998) 

For  a  PMMA  sphere  at  largest  resonance 

(to  =  1.73) 

a{m)  =  AOllF(Hz) 

R„f=  1  m 
1/1=5.63 


TS  =  20  log 
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TS  (dB) 
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Lamb  waves  on  an  acrylic  plate. 

Phase  velocities  for  2mm  thick  acrylic  plate 
without  fluid-loading 


Characteristic  equation  for  symmetric  waves  on  a  plate  immersed  in  water. 
Vacuum  term  cosh(aA/2)sinh(p/i/2)-4^^aPcosh(pA/2)sinh(a/i/2)J(- 


Fluid  -  loading 


„2  e2 


-^^^inh(pA/2)sinh(a/i/2) 


=  0 


Two  properties  of  acrylic  which  significantly  affect 
the  fluid-loaded  Lamb  wave  spectrum: 


1)  Cr<Cw 

2)  Pe  ^  Pw 
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For  Acrylic: 


P.  ^  1.00 
Pe  1.19 


0.84 


For  Stainless  Steel:  =  4^  =  0.13 

Pe  7.57 

For  large  density  ratios,  fluid-loading  is  no 
longer  a  small  perturbation  to  Lamb  wave 
spectrum. 


As  pw/pE  increases,  the  Lamb  wave  spectrum  transitions  to 
a  plate  with  mixed  boundary  conditions. 

S.l.  Rokhlin  et.al.,  J.  Acoust.  Am.  85(3),  March  1989 
A.  Freedman,  J.  Acoust.  Am.  99(6),  June  1996 
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Damping  parameters  for  an  empty  fluid-loaded  acrylic  shell 


Summary 


•  PMMA  spheres  were  observed  to  give  a  backscattering 
enhancement  due  to  tunneling  to  a  subsonic  Rayleigh 
wave. 

•  We  have  modeled  the  resulting  Rayleigh  wave  resonance 
scattering  by  extending  the  previous  ray-tunneling 
models  of  subsonic  ao  waves  on  steel  shells  [Zhang,  Sun, 
&  Marston  JASA  (1992)]. 

•  Radiation  damping  and  intrinsic  material  absorption  are 
similar  in  magnitude  and  are  additive. 

Computational  study  for  shells: 

•  Because  the  Rayleigh  wave  velocity  is  subsonic,  there 
are  two  branches  associated  with  the  lowest  symmetric 
mode,  the  So  and  Sq.  waves  which  give  rise  to  associated 
backscattering  enhancements. 

•  Although  the  intrinsic  material  absorption  may  suppress 
enhancements  associated  with  waves  which  travel 
around  the  shell,  it  may  be  possible  to  detect 
enhancements  associated  with  a  backwards  wave  on  the 
shell. 
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Noninvasive  Identification  of  CW  Agents 
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Remote  Classincation  of  Fluids 
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Firefighters,  HAZMAT 


Noninvasive  Identification 
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situation  (e.g.,  chemical 
vapor,  radioactive,  etc.,) 


Nonlinear  Sound  Propagation  in  Air 
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Remote  Excitation  of  Vibration 
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Remote  Classification  of  Fluids 
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Beamwidth 


Directivity  of  Parametric  Array 
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Degrees 
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Laser  Vibrometer 


Instrument 


/ — \ 
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Laser  Vibrometer 


Ultrasonic  Remote  Assay  of  Munitions 
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containers 


Liquid-filled  Steel  Containers 
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•  interior  fluid  attenuation 

•  shell  diameter 

•  shell  thickness  or  material  composition 


Calculated  Phase  Velocity  and  Radiation  Damping 

for  a„  Lamb  Wave 


Radiation  Damping  (Np) 
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Frequency  (kHz) 
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Resonance  Spectrum  of  Two  Liquid-filled 

Steel  Containers 
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Frequency  (kHz) 
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It  is  possible  to  monitor  liquid  level  remotely 

It  is  also  possible  to  determine  wall  integrity 

This  remote  monitoring  approaeh  has  many  applications 


Acoustic  Gas  Resonators  for  Measurements 
of  Thermophysical  Properties  and 
Thermometry 


James  B.  Mehl 


Resonance  Meeting,  1999 
University  of  Mississippi 


University  of  Delaware 
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u  +  iv\  vs.  / 


/o 
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Interesting  parameters:  F  =  /o  + 


u  +  iv  =  — — —  +  background 

jz  _  p2 

Insert  extra  factor  of  /  in  numerator: 


/  =  F  +  (/-F)  =  F  +  A/ 


Then 


iA/F  iAfAf  ,  , 

“  +  +  (j  +  background 

iA'f  (l  + 

=  - ^ — A~f~^  background 


iA'f 


+  modified  background 


r  -F2 

Similar  relations  will  handle,  e.g. 


A{f)  =  A(F)  +  ^Af  +  ... 
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Frequency  (HZ) 

FIG.  4.  Top:  In-phase  (solid  curve)  and  quadrature  (dashed  curve)  vol¬ 
tages  from  the  detector  as  a  function  of  frequency  near  the  (0,2)  resonance 
in  argon  at  0.4032  MPa  and  296.309  K.  Bottom:  Measured  voltages  minus 
calculated  voltages  (Eq.  (75)  ]  with  the  fitted  parameters /02  =  3613.9970 
Hz,  go2  =  0.4275  Hz,  A  =  0.7535  -  0.5831/,  B  =  0.0048  4-  0.0068/,  and 
C=  0.000 i  +  0.0001/.  Data  were  taken  at  intervals  of  0.07  Hz,  first  from 
the  lowest  to  the  highest  frequency  and  then  back  to  the  lowest  frequency. 
Separate  deviation  plots  for  the  upward  and  downward  frequency  sweeps 
show  the  effects  of  drifts  in  the  lock-in  amplifier. 
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FIG.  5.  Top:  Ih-phase  (solid  curve)  and  quadrature  (dashed  curve)  vol¬ 
tages  from  the  detector  as  a  function  of  frequency  near  the  (1,2)  resonance 
in  argon  at  0,4032  MPa  and  296.309  K.  Middle:  Measured  voltages  minus 
two-resonance  trial  function.  Note  that  the  deviations  are  systematic  al¬ 
though  the  trial  function  has  twelve  parameters  [eight  parameters  specify 
resonances  at  4777.63,  and  4779.38  Hz,  and  four  specify  the  constant  and 
linear  background  terms  in  Eq.  (75)].  Bottom:  Measured  voltages  minus 
fitted  function.  The  fitted  function  has  sixteen  parameters  [twelve  param¬ 
eters  specify  resonances  at  4777.693, 4777.903,  and  4779.351  Hz  with  half¬ 
widths  of  0.550, 0.546,  and  0.555  Hz.  The  remaining  four  parameters  spe¬ 
cify  .5  and  C  in  Eq.  (75)].  jj.. 
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Figure  5.  Cross-section  of  resonator  and  pressure  vessel.  The 
transducer  assemblies  are  indicated  by  “T,”  and  the  locations  of 
the  capsule  thermometers  are  indicated  by  “PRT.” 
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Acoustic  eigenvalues  of  rigid,  insulating  sphere 


A  gas  filled,  acoustic  resonator  with  rigid,  thermally  insulating  walls,  has  acoustic  resonances  at 
the  frequencies 

fin  -  ^  (0.577  kHz)2:M, 

where  the  numerical  value  is  for  argon  at  room  temperature  (c  «  322m/s)  in  the  8.89  cm-SS 
spherical  resonator.  The  acoustic  pressure  in  the  sphere  is  proportional  to 

3i{ztnr/a)Y(,n{0,(t>). 

The  modes  are  +  l)-fold  degenerate  for  perfect  spherical  geometry. 
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Theory  of  spherical  acoustic  resonator 

Temperature,  pressure,  particle  velocity  fields 
PDEs  (Kirchhoff  1868) 

Navier  Stokes:  longitudinal 
Navier  Stokes:  transverse 
Fourier:  Heat  flow 
Conservation  of  mass 
Conservation  of  energy 

4th  order  PDF  in  temperature 
Prediction:  two  wave  modes 
Ap,  AT,  particle  velocity,  oc 
acoustic  mode  {kp  ^  u/uj) 
thermal  mode  {kt  (1  -  ^)/St) 

Boundary  conditions  at  gas-shell  boundary: 

velocity  of  gas  =  velocity  of  shell 
heat  flow  continuous 

temperature  continuous  (hydrodynamic  region) 

Shell  motion: 

Classical  elastodynamic  theory  exactly  soluble 
for  isotropic,  isolated,  spherical  shell 

Complete  problem: 

Exactly  soluble  for  spherical  geometry 
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Acoustic  experiment 

1.  Measure  resonance  frequencies  finm  and  half-widths  gimn 

2.  Correct  for 


(a)  thermal  boundary  layer 

(b)  viscous  boundary  layer  (non-radial  modes) 

(c)  shell  motion 

(d)  imperfect  geometry 


^fv  +  i9v  _  ^  -I  Sy  l(l  +  1) 

/  “  2  a  I- 1(1 +  1)1  zl 
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Boundary  shape  perturbation: 


r  =  a  [1  —  eT{6,  (j))] 


e  =  deformation  parameter. 


For  deformations  that  conserve  the  volume; 


[TJon 


0(6^) 


It 


=  0(6) 


For  nonradial  acoustic  and  electromagnetic  modes: 

Average  over  2^+1  components  of  a  multiplet. 


/ 


geom  =  O(e^) 

tn 
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Shell  admittance  function 


Shell  admittance  function  for  the  stainless-steel  resonator  with  a  =  8.89  cm  and  h/a  =  1.215 
(a  =  inner  radius,  b  =  outer  radius).  The  curves  are  labeled  with  the  mode  index  £,  for  acoustic 
pressures  proportional  to  Pe{cos9).  The  fractional  frequency  perturbation  due  to  shell  motion  is 

V  /  /  ~  1  -  +  l)/zj^  ■ 

Here  =  ujjcsh^i,  with  Csh,i  the  longitudinal  speed  of  sound  in  the  shell. 
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Frequency (Hz) 

FIG.  1 1 .  Response  of  the  shell  to  radially  symmetric  excitation  as  a  function 
of  frequency.  The  points  are  the  average  slopes  of  the  curves  in  Fig.  9.  The 
curve  is  calculated  for  an  isotropic  seaniless  shell  using  the  theory  of  elasti¬ 
city  and  the  elastic  constants  tabulated  for  aluminum  (see  Table  V).  The 
idealized  shell  has  a  breathing  resonance  near  20.2  kHz. 
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Frequency  (Hz) 

FIG.  12.  Elastic  response  of  the  shell  to  excitation  with  symmetry  of 
Yim  as  a  function  of  frequency.  The  points  are  the  average  slopes  of 
curves  such  as  those  shown  in  Fig.  10.  The  curve  is  calculated  for  an  isotrop¬ 
ic,  seamless  shell  using  the  theory  of  elasticity  and  the  elastic  constants  tabu¬ 
lated  for  aluminum.  The  idealized  shell  has  resonances  at  0  and  24.5 
kHz. 
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Sampled  Continuous  Wave 
measurement  on  single  crystals  of 
YBCO  and  BSCCO. 


Debashis  Dasgupta 
Jeff  Feller 

Carsten  Hucho 
Moises  Levy 
Bimal  K  Sarma 


Resonance  Meeting,  1999 
Univ.  of  Mississippi 


Graduate  School,  UW-Milwaukee 
Physics  Department,  UW-Milwaukee 
Office  of  Naval  Research 
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Small  size  of  HTSC  single  crystals, 
(typically  ~  1000  X  800  X  200  microns) 


Objective; 

Develop  an  ultrasonic  technique  for 
investigating  properties  of  these  small 
HTSC  crystals. 

Sampled  CW  Technique 
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OUTLINE 


Type-II  Superconductors 

Essentials  (Magnetic  phase  diagram) 
Flux  lines  and  Pinning  sites 
Description  of  Sampled  C  W  Setup 
Experimental  Results  On  YBCO 
Untwinned  Crystal 
Twinned  Crystal 
Susceptibility  Results  (Twinned) 
Experimental  Results  on  BSCCO 
Conclusions 
Future 
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C2 


A  schematic  of  complex  magnetic  phase  diagram  of  a  Type-II 
superconductor.  Some  of  the  regimes  shown  are  magnified  for  clarity. 
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3)  Magnus  Force 


S 

Moving  vortex 
(Lab  frame) 


Static  (spinning) 
(vortex  frame) 


medium 


velocity  - 


Pinning  sources  for  flux  lines 


Point  defects 

•  Oxygen  Vacancies 

•  Electron  Irradiation 
induced  point  defects 


Line  Defects 

•  Heavy  ion  irradiation 
induced  columnar  tracks 


Planar  Defects 

•  Twin  Boundaries 

•  Layered  Structure 
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Continuous  Wave  Technique 


>•  Transducer  < 


Sample  length  x 


•  Setup  standing  wave 
resonances 


Q 


m 


CO 


m 


Aco 


Lock  in  to  mechanical  resonance 

Change  external  parameters  (Temperature,  Magnetic  Field  ... ) 


Af  ►  Av  ^  Aa 
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Sampled  CW  Experimental  setup:  transducer  +  crystal 

B 


X  »  thickness  of  the  sample  assures  no  displacement  gradient  for  the  ions 

sample  thickness:  ~50|xm  obtained  through  polishing  (optical  microscope  for  thickness) 

sound  velocity  c:~  5000  m/s,  frequency:  5MHz  X~  I  mm  »  sample  thickness 


Temperature  (K) 


IIq2  and  melting  curves  for  a  twinned  YBCO  crystal. 

( H  //  c-axis) 


in  dB/fisec 


Conclusions; 


Sampled  CW; 

•  “The  ultrasonic  technique” 

•  Minimum  external  perturbations 

thermodynamic  measurements 

•  Bulk  elastic  properties  of  flux  lines 

•  melting  of  flux  lines  in  YBCO  and 
BSCCO 

•  soft  lattice  at  low  fields 

•  rigid  lattice  at  higher  fields 

•  possibly  soft-rigid  contour  follows 
the  melting  curve  with  0.3T  as  a 
turnaround  field 

•  Sensitive  to  the  “depinning”  transition 

•  change  in  activation  energy  around 
0.3T  for  YBCO 

•  Basic  electron-phonon  interaction 

SC  transition  at  11^^ 


Susceptibility; 

•  Peak  effect  for  YBCO 

•  Re-entrant  nature  of  the  peak  effect  has  been 
observed 
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